Spin injection into a metal from a topological insulator 
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We study a junction of a topological insulator with a thin two-dimensional (2D) non-magnetic 
or partially polarized ferromagnetic metallic film deposited on a 3D insulator. We show that such 
a junction leads to a finite spin current injection into the film whose magnitude can be controlled 
by tuning a voltage V applied across the junction. For ferromagnetic films, the direction of the 
component of the spin current along the film magnetization can also be tuned by tuning the barrier 
potential Vb at the junction. We point out the role of the chiral spin- momentum locking of the Dirac 
electrons behind this phenomenon and suggest experiments to test our theory. 

PACS numbers: 73.20.-r,73.40.-c,71.10.Pm 



Topological insulators (TI), a class of three- 
dimensional (3D) insulators with strong spin-orbit 
coupling, are known to possess gapless Dirac-like quasi- 
particles on their surfaces whose existence originates 
from the special topological properties of their bulk 
bands [l|49|. Such insulators, which are essentially 3D 
generalizations of their 2D counterparts which exhibit 
the quantum spin-Hall effect, have attracted a lot 
of theoretical and experimental attention in recent 
years. These TIs can be classified as strong or weak 
depending on their sensitivity to time reversal symmetric 
perturbations. The surfaces of the strong TIs have an 
odd number of Dirac cones; the number and positions 
of these cones depend on the nature of the surface 
concerned [l|, |2|, Ij] . The odd number of the Dirac cones 
ensures that any surface impurity which conserves time 
reversal symmetry does not destroy the low-energy 
Dirac properties of the quasiparticles on the surface. For 
several compounds such as Bi2Te3 and Bi2Se3, specific 
surfaces have been found with a single Dirac cone near 
the r point of the 2D surface Brillouin zone [2|, [Zl |9j . 

A Dirac cone on the surface of a TI is described by the 
Hamiltonian 




FIG. 1: (Color online) Schematic picture of the proposed junc- 
tion. The two interfaces I and II of the topological insula- 
tor have low-energy quasiparticles with Dirac-like properties 
while the non-magnetic/ferromagnetic metallic thin film III 
deposited atop the insulator has conventional electrons obey- 
ing the Schrodinger equation. See text for details. 



tance of ferromagnet-ferromagnet junctions [I4I, realiza- 
tion of a magnetic switch in junctions of these materials 
[13| , spin textures with chiral properties [9|, control of 
spin transport and polarization using gate voltages and 
"10, 



Hnivp] = 



dkidk 



^ ipi (hvpfi -a xk- ^I) tpj:, (1) 



where (?(/) denotes the Pauli (identity) matrices in spin 
space, n denotes the unit vector normal to the TI sur- 
face which hosts electrons with momentum components 
ki and kj, vp is the Fermi velocity, /i is the chemical 
potential, and ip = (V't, V'i)'^ is the annihilation oper- 
ator for the Dirac spinor [10[. In this notation, t {I) 
denotes components of the quasiparticle spin along (op- 
posite to) z. Recently, several novel features of these 
surface Dirac electrons have been studied. These include 
the existence of Majorana fermions in the presence of a 
magnet-superconductor interface on the surface 10|, [ll [ , 
generation of a state resembling apa;-|-ipy-wave supercon- 
ductor but with time reversal symmetry via proximity to 
a s-wave superconductor 'ICr], anomalous magnetoresis- 



electric fields 1 14 Il5t , realization of a Lifshitz transition 
in thin TI films [l6| , and spin-polarized STM spectra [l7[ . 
However, junctions of such TIs with conventional metals 
and ferromagnets have not been studied so far. 

In this Letter, we study the transport properties of 
a junction of a TI with a conventional insulator with a 
non-magnetic or partially polarized ferromagnetic metal- 
lic film deposited on it as schematically shown in Fig. 
[T] The regions I and II in Fig. [T] refer to the two sur- 
faces of the TI which host chiral Dirac quasiparticles, 
while region HI consists of conventional electrons obey- 
ing the Schrodinger equation. We show that due to the 
chiral spin-momentum locking of the Dirac electrons on 
the surfaces of the TI, the transport through such a junc- 
tion can lead to a finite spin current in the film without 
any skew scattering which is necessary in conventional 
spin-Hall materials to generate such currents. The mag- 
nitude of the spin current generated can be controlled by 
a voltage V applied across the junction without the in- 



volvement of any external magnetic field. We also study 
the dependence of this spin current on the barrier po- 
tential Vq at the junction and show that it displays an 
oscillatory feature with a decaying envelope as a function 
of Vq. This behavior draws from both the chiral nature 
of the Dirac quasiparticles in regions I and II which leads 
to the oscillatory nature of the spin current [18| , and the 
presence of the conventional Schrodinger electrons in re- 
gion III which leads to an exponential decay of the spin 
current with increasing Vq- Finally, we demonstrate that 
for ferromagnetic films, the direction of the component of 
the spin current along the film magnetization {J^) can be 
controlled by tuning Vq. To the best of our knowledge, 
the generation of spin current using the chirality of the 
Dirac quasiparticles on the surface of the TIs whose di- 
rection and magnitude can both be controlled electrically 
has not been proposed before; we therefore expect our 
proposal to generate significant interest in the field of 
spintronics. 

We begin with the analysis of the junction in Fig. [T] 
when Vb = 0. In region I, the Hamiltonian for the Dirac 
quasiparticles is given by Eq. ([l} with h = z. The wave 
function for these quasiparticles with a transverse mo- 
menta kx and energy e — eV moving along ±y can be ob- 
tained by solving the Dirac equation Hz[vi\'ipi = e^'i and 
is given by V'± = [1, -i exp(±j6'fc)]^ eyiY>[i{kxX±k yy)]/ \/2 , 
where 9^ = a,rccos[hvikx/{e~\-^)], e — —ii+hvi^/k'^ + ky, 

and vi is the Fermi velocity of the quasiparticles in re- 
gion I. The wave function of the Dirac quasiparticles in 
region I can thus be written as 



'ipi = ip++ rip- , 



(2) 



where r is the refiection amplitude. We note that for 
any incident angle 0, ip+ and -0- have {ip+\ax\ip+) = 
sin(6'fc) = —('(/'_ I tJa; IV'-)- Thus the refiected Dirac quasi- 
particle has the opposite spin orientation along x com- 
pared to its incident counterpart. This phenomenon is 
reminiscent of Andreev reflection from a superconducting 
junction where the charge and the transverse momenta 
of the refiected quasiparticle change sign. In contrast, 
for the junctions considered here, the transverse in-plane 
component of the quasiparticle spin changes sign upon 
reflection. In what follows, we shall show that this spin 
reversal is at the heart of the generation of a flnite spin 
current in region III. 

In region II, fi = y, and the wave function of the Dirac 
quasiparticles moving along —z with transverse momenta 
kx and energy e can be obtained by solving Hy[v2]ijj2 — 
£■02 and is given by 



tpn = tl'02, ip2 = [uk,Vk 



iT i{-k^z+k^x) 



2, (3) 



cos 



c): 4>k 



where Uk[vk] = y/T^ 

aTCCos[hv2kx / {e + n)], e ~ -fi + hv2y/^+k^, V2 = P'^vi 
is the Fermi velocity, and ii denotes the transmission 



probability of the Dirac quasiparticles in region II. In 
the rest of this work, we shall choose j3 = \Jv2lv\ < 1. 

In the metallic film (region III), the Hamiltonian 
for the electrons can be written as Hjjj — fi^{kx + 
ky)/{2m) — fJ. — Aaz, where /i and m are the chemical 
potential and mass of the electrons in the film, and A is 
proportional to the magnetization of the electrons. For 
a non-magnetic film A = 0, while for a fully polarized 
ferromagnetic film A — > oo. In what follows, we first 
consider a non-magnetic film for which ^ = 0. The wave 
function of the electrons in region III is then given by 



t 



III 



[t2,t3f e*(fe-^+fe;!') 



(4) 



where e = —fi + h?{kx + ky)/{2m), and ^2 and is denote 
the transmission amplitudes of spin-up and spin-down 
electrons in region III. 

The boundary condition on these wave functions in- 
volves continuity of current through the junction which 
yields 



vitplaxTpi - V2ip\iax'ip 



II 



Im{ipljjdyipiii), (5) 



where it is understood that all fields are evaluated at 
the junction line y = (for ipj and ipm) and z — Q 
(for V'//)- We note that the unusual boundary condition 
(Eq. [5]) in which a current without derivatives in regions 
I and II (Dirac equation) has to be matched with a cur- 
rent involving a first derivative in region III (Schrodinger 
equation) is generic for any junction involving a TI and 
a conventional material. (The situation here is different 
from a junction of ordinary materials with spin-orbit cou- 
pling where the current involves a first derivative on both 
sides [19|). Below, we present a general solution to this 
problem with (discussed in Eq. p^ below) and without 
a barrier potential 20] . 

In the absence of any barriers at the junction, the gen- 
eral solution of Eq. ([5|) is given by two linear conditions 



on the wave functions J21| . 



-0/// = c(V'/ + p-ipii), —dyipin 
m 



IVlCTx 



l^i - Pt 



II 



(6) 

where c is an arbitrary real constant; we will set c — 1 for 
simplicity. Note that the metal/ferromagnet decouples 
from the TIs for c — > or oo. Substituting Eqs. (|2]|4|) in 
Eq.[6l we obtain the following relations between r, ti, t2, 
and ta, 

l + r+{-)l3ukti = t2iat3), 
£'»'=+ re-'«'= + (-)i/3i;fcti = itz{iat2), (7) 

where a — hk' / {mvi). Solving for r, ii,2,3 from Eq. ([7]), 
we get 



r - Af/V, ti - 2(1 - a") sm{0k)/{(3V), 



(8) 



4 4 

t2 = ^sin(6'fc)(ufc -l-aufc), t^ ^ —sm{9k){vk + auk), 
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FIG. 2: (Color online) Plot of J^ (top panel) and I^ (bottom 
panel) in a non-magnetic metallic film as a function of V for 
/3 = 1. The insets show their /3 dependence for eV jEt^ = 3. 



where A^ = — zMfc[exp(i^fe)(l -f a^) — 2ia\ — Wfc[(l -l-a^) + 
2iaexp(j0i;)], and V = i[Mjt(l + a^) + 2awfc] exp(— 10^) + 
Wfc(l -f a^) +2aMfc. 

Eq. ([HI) predicts a net spin current along x in the 
metallic region: J^ = (St^i/2) X]fe^(V'///|at^a;|'0///) = 
{hvi/2) Y^j, a[t2t^ -f h.c). This can be written as 



, _8Jo 

TT 



E/{2Ea) 



dx 



E/(2Eq) 



t[(l + a^)sin(0fc) + a]sin^(6>fc) 

(9) 

where Jg = fivikg, x — k^/ko, kg = mvi/h, E = eV + /i, 
and Eq — hvikQ/2. A plot of Jx/Jo as a function of the 
applied voltage eV/Eo, shown in the bottom panel of Fig. 
[51 confirms that the net spin current is finite and its am- 
plitude depends on V. The inset shows the dependence of 
the spin current on /3 for a fixed V and confirms the pres- 
ence of a finite Jx for the entire range of < /3 < 1. The 
charge current is given by Ic = evi J2k Q^(l*2p + l^sP) 
which can be written as 



4io 

TT 



E/i2Eo) .^ 

dx 

'E/i2Eo) 



a^ + 2Qsin((/)fc)]sin^(6>fc) 
IPP 



(10) 

where /q = eviko. The top panel of Fig. 2 shows the 
variation of Iq/Iq with the applied voltage V, while the 
inset depicts its variation with /3. We note that the 
charge current displays a qualitatively similar behavior 
as the spin current along x. The net spin current along 
y and z vanishes. For J^ this can be seen by noting that 
Uk — > Vk under k^ -^ —kx- Consequently t2 — > t^ un- 
der this transformation which leads to a zero net value 
for Jz ~ J2k o^d^sP — l^sP)- Also, using the fact that 
t2 and ^3 are both purely imaginary (Eq. (|8l)), it can be 
easily shown that Jy = [22| . 

Next, we address the behavior of the spin and charge 



FIG. 3: (Color online) Plot of Jz versus eV when n = A (red 
dashed line) and /i — —A (blue solid line) for A/Eo = 3 (top 
panel) and as a function of A (bottom panel) for P — 1 and 
eV/Eo — 3 and n — A (red dashed line) and fj. — ~A (blue 
solid line). 



currents for a ferromagnetic film where A ^ 0. In this 
case, the wave function in region III is given by 



V'. 



FM 
III 



[t2e<'\t3e<''y] 



T ik^x 



(11) 



(1)[(2)] 



where 4'^^^^^' = [2m{e + ^i + [^]A) / h^ - kl\' . The 
boundary condition on the wave function is given by 
Eq. ((HI and yields 

^2(3] == -4ie^'''=sin(6ife)(ufe[a2Ufe]+aif/cK])/2?i,(12) 

where ai(2) = hky ' [ky ')/{mvi), and Vi — Pi — 
ip2exp{i9k)j where pi — Ufe(l -I- aia2) + 2vka2, and 
P2 — Vk{l + aia2) + 2ukOt2- Note that for ^ = 0, 
cti = a2 = a; in this limit Eq. ([12]) matches with Eq. ([8|). 
Using this wave function, it is straightforward to com- 
pute the expression for Ic following the method outlined 
earlier which shows qualitatively similar behavior as in 
metallic films. 

The key difference between the ferromagnetic and the 
non-magnetic films which we now focus on is that the 
former films allow a non-zero Jz- This is easily seen 
from Eq. (fT2|) by noting that t2{kx) ^ H{—kx) due to 
the difference of velocities of the up and the down spin 
quasiparticles; this leads to a finite J^ in region III given 
by 



J. = 



J j-E/{2Eo) 
•''0 I j„^ /„ u |2 „ I, |2\ 



27r 



dx{a,\t2\'-a2\t3n. (13) 



-E/{2Eo) 



In what follows, we set /i = A or /i = — A so that the 
spin-up or spin-down Fermi surface is aligned with the 
Fermi surface in region I. Increasing A therefore pushes 



the other Fermi surface away from the Fermi surface of re- 
gion I and hence increases the effective spin polarization 
of the film in region III. The behavior of J^, computed 
by substituting Eq. p^ in Eq. (IT^ . is shown in Fig. [31 
The top panel of the figure shows the variation of Jz as a 
function of cV/Eq for A/Eq ~ 3, while the bottom panel 
shows the dependence of Jz on A for fixed eV/Eo ~ 3 
and (3 — 1. These plots demonstrate the presence of a 
finite Jz for a large range of V; the sign of Jz depends 
on the magnetization of the film while its magnitude can 
be controlled by the applied voltage V. 

Next, we consider the effect of a finite barrier potential 
Vq applied over a region d at the junction. In what follows 
we will consider the thin-barrier limit where Vq ~^ oo and 
d — > keeping x = Vod/ihvi) finite. In this limit, the 
boundary condition to be imposed on the wave functions 
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4^111 = e-^^"==V^/+/?e^^"==/^V//, (14) 
-dy^iH -2x4^111 = icjAe-'^^-i,i-^e^^-lf"4'u)- 



Substituting Eqs. ([2]), ([3]), and dlT]) in Eq. dUl), we again 
obtain a set of four equations for r and ti.2.3 which reads 



t 
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FIG. 4: (Color online) Top left (right) panel: Variation of J^, 
(/c) with X in a metallic film for eV jE^, = 3. Bottom left 
(right) panel: Variation of J^ with x in a ferromagnetic film 
for /i = A = 3So (ai = -yl = 3So) for eVjEo = 4 (red solid 
line), 3 (black dotted line) and 1 (blue dashed line). For all 
the plots /3 = 1. 



To compute the spin and charge currents, we sub- 
stitute the values of ^2(3) from Eq. (fTS]) in the expres- 
sions Jz = ^wi X]fc (Q^i|'^2p — a2|i3p) for ferromagnetic 



(ai[m2]4-ra*[m;] + [-]/3ii&i[62]), films, and I^ = (^v\Y.k ai(l*2p + 1*3^) and J^ 



("1(2] + 2ix)^3[2] = -;^(ai[-*a2] + ra\\-ia*^ - /3ii6i[fo2]) 

where a^ — cos(x) — sin(x) exp(i0fe), a^ — sin(x) + 
cos(x)exp(i0fe), and 6i[2] = Ufe[wfe] cos(x/;S2) + 
JWfe[ufe] sin(x//3'^)- The solution of these equations yields 

i2[3] = [-2(5i[2] + &2[i](ai[2] + 2ix))(a\a2 - alai)]/V2, 

(15) 

where r'2 = [-26i(2x -I- JQi) + &2{-l + ("i - 2ix){a2 - 

2zx)}]mi + [6i(l + (ai+2zx)(a2 + 2zx)) + 262(a2 + 2ix)]a2 
for magnetic films. The corresponding expressions for the 
non-magnetic films can be obtained by putting ai = 02 
in Eq. ^. We note that Eq. ([TS]) reproduces Eq. ^ 
for X = 0. 

From Eq. (jlSp , we find that the barrier potential x en- 
ters the transmission amplitudes ^2(3] both through the 
cos(x) and sin(x) factors in 01(2) and 61(2) leading to an 
oscillatory x dependence of ^2(3) and through the appear- 
ance of X in the numerator and denominator of Eq. (|15p , 
which, in the limit of large x, leads to a decay of i2(3) with 
increasing x- The former behavior arises from the Dirac 
nature of the electrons in regions / and //, while the 
latter is a consequence of the conventional Schrodinger 
nature of the electrons in region III. Consequently, we 
expect ^2(3) to have an oscillatory dependence on x along 
with a decaying envelope. We note that such a behav- 
ior is different from what is found in analogous junctions 
involving solely Dirac or solely conventional materials. 



^"1 Sfc 0^(^3*2 + h.c.) for non-magnetic films. The re- 
sulting dependence of Jz, Ic and J^ on x is shown in 
Fig. 0] for (3 = 1. The top panel shows the behavior of 
Jj, and Ic for non- magnetic films for cV/Eq — 3. We 
find that both J^ and Ic display small oscillatory fea- 
tures with an overall monotonic decay as a function of 
X- We have found qualitatively similar behavior of Jx 
and Ic in magnetic films. In contrast, Jz for magnetic 
films, shown in the left (right) bottom panels of Fig. 0] 
for li = A {ji = —A) for several representative values of 
eV/ E(j displays a non-monotonic behavior with increas- 
ing X- In particular, these plots demonstrate that the 
sign of the spin currents gets reversed with increasing x 
which allows for the possibility of tuning the sign of Jz 
by tuning Vq. This leads to electrical control of both the 
magnitude and the direction of Jz in magnetic films via 
tuning externally applied voltages V and Vq. 

To experimentally verify our theory, we propose mea- 
suring the current in region III for magnetic films using a 
ferromagnetic contact. When the direction of magnetiza- 
tion of the contact is along z, it will measure only the cur- 
rent due to the spin-up electrons: I^ = ewi ^j, ai|i2p — 
[Ic -t- eJz/(2/i)]/2. Similarly a contact with magnetiza- 
tion along —z will record /j. = [Ic — eJz/{'2,h)]/2. Thus 
the difference between these two currents will provide a 
measure of the spin current through the film as a func- 
tion of the applied voltage V and the potential barrier 
Vq: Jz = 2/i(/-f- — Ii)/e. Thus such an experiment can 
verify the predicted dependence of Jz on V and Vb [2l| . 

In conclusion, we have studied a junction of a TI with 



a thin metallic/partially polarized ferromagnetic film de- 
posited over an ordinary insulator. For ferromagnetic 
films, we have shown that such a junction can be used 
to generate a finite spin current along z whose magni- 
tude and direction can both be controlled by externally 
applied voltages without the presence of any external 
magnetic field. Our work shows that the chiral spin- 
momentum locking of the Dirac quasiparticles on the 
surfaces of the TI is at the heart of this phenomenon. 
Finally, we have suggested a simple experiment to test 
our theory. 
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APPENDIX A 

In this Section, we provide some supplemental material 
to the main text related to the derivation of the boundary 
condition Eq. [TJ] and numerical estimate of the measured 
current in suggested experiment. 



Current conserving boundary conditions at a 
junction 

In this section, we will find the general time reversal 
invariant boundary condition which satisfies the current 
conservation relation at the junction discussed in our pa- 
per. We begin with the Hamiltonians in the three regions. 
Region I of the topological insulator (TI) is defined by 
the half-plane z = Q and j/ < 0, and has the Dirac Hamil- 
tonian 



Hi = ihvi [~axdy + cFydx]- 



(16) 



Region II of the TI, given by the half-plane y = and 
z < 0, has the Hamiltonian 



H 



II 



ihv2 [-(T^dx + a^d;, 



(17) 



Region III of the non-magnetic metal/ferromagnet thin 
film is defined by the half-plane z = and y > 0, and 
has the Schrodinger Hamiltonian 



H 



III 



2m 



{dl + dl) - M - Aa.. (18) 



The time evolution equations ihdtpa/dt = Haipa (a = 
/, //, ///) are invariant under time reversal (i — > —t and 



complex conjugation of all numbers) if 4'i ~^ '^y'^t-'t ^^d 
A = 0. 

The junction of the three regions is given by the line 
y = z = Q. The conservation of the total current coming 
into the junction from the three regions is given by Eq. 
(5) of our paper, namely, 



',t. 



',t 



= —Im{ii\jjdy'il)iH)y=Q+. (19) 

In order to satisfy this equation, let us assume linear 
relations between the wave functions at the junction of 
the form 

ilplll)y=0+ == Ai{'lJji)y^o^ + A2{ll^Il)z=0-, 

— idyil;iii)y=o+ 



= ia-x [A3{iJJi)y=o- + A4('0//)2=o-], 

(20) 

where the A^ are four parameters. The relations in 
Eq. (PO)) will be time reversal invariant if the Ai are 
real. We can now check that Eq. (|19p will be satisfied 
if A1A3 = vi, A2A4 = ~V2 and A1A4 + A2A3 = 0. This 
implies that the Ai can be written in terms of a single 
real parameter c as Ai = c, A2 = c/3, A3 = vi/c and 
A4 — —viP/c, where /3 — \/ V2lv\\ this gives 



('0///)y=o+ 

-(9yV'///)y=0+ 



c[(7/'/)y=0- +^(V'//)z=0-], 
—^-^ [(V'/)y=0- - /5('0//).z=O~], 



(21) 



which is Eq. (6) of our paper. In the limits c -^ (or 
00), we obtain ■(/;/// = (or dyipni =0); in either case, 
the current into the junction from region III vanishes, so 
that the metal/ferromagnet gets decoupled from the two 
TI regions. The value of c in a given system will depend 
on its microscopic details such as an underlying lattice 
model. For the numerical calculations in our paper, we 
have simply set c = 1. 

The above analysis assumed that there is no barrier 
present at the junction. A realistic system may be ex- 
pected to have some potential barriers present at the 
junction in all three regions. Let us assume thin bar- 
riers of the form V{y) — Vq for ~d < y < in region I, 
V(z) — Vq for — d < z < in region II, and V{y) — Vq for 
< y < d in region III. For simplicity, we have assumed 
the barrier width (d) and height (Vq) to be the same in 
all three regions; we will eventually be interested in the S- 
function limit d — > and Vo ~^ c» keeping dVQ/{hvi) = x 



constant. In Ref. ? , it has been shown that 5- function 
barrier in a Dirac Hamiltonian produces a discontinuity 
in the wave function of the form 

(V/),=o- = e-*x-^ {^Ji)v=~d, 

{.^ii)z=o- = e^>^("^/''^)-^ (V'//).=-d. (22) 

A (5-function barrier in a Schrodinger Hamiltonian 
produces no discontinuity in the wave function (i.e., 
('0///)i/=d = {4'iii)y=o+), but there is a discontinuity in 
the first derivative of the form 

— [{dy%PlIl)y = d - {dyi>IIl)y^O+] = 2X^1 ( V'/// ) ^ = d ■ 

(23) 
Substituting Eqs. (1221123^ in Eq. ([21]), and setting c = 1, 
we obtain 



iiJlIl)y=d 



— p-^X<^L 



(V'/)y=-d + /3e*^'^/^'^""'(V'//).=-d, 



— {dy'4>ui)y=d - 2xwi('(/'///)a=d 

= max[e-*'^"-(V'7)j,=-d-/3e^^x/^'^"-(^//).=-d]. (24) 
In the limit d — >■ 0, this gives Eq. (14) of our paper. 

Numerical estimates of measured currents 

In this section we provide estimates for I^ and I^ which 
are to be measured in the proposed experiment. For a 
typical TI surface, the group velocity of the Dirac elec- 
trons turns out be vi ^ lO^m/s. This give us an estimate 
of /co — rnvi/h ~ 8.63 x 10~^m~^. Using this, one finds 
/o = evikQ ~ 13.83mA and £^0 = mvl/2 ~ 2.89eV. Since 
I^ and I^ can be written as 



4/„ rE/{-2Ea) 



'T(i) 



TT J-E/{2Eo) 



dxailial^ (aalial'), (25) 



we find numerical values of /-f- = 0.113mA and I^ = 
0.0906mA for x = and fi = -A = SEq. This indi- 
cates that the visibility V defined by 



V = 



h - h 



h + h 



(26) 



is close to 0.1 which means that such current measure- 
ments are well within the reach of current experimental 
standards. For finite and large barrier strength x = 5, 
the corresponding numbers are /-|- — 0.0068mA and 
4 = 0.012mA which leads to F ~ 0.28. 



